Finite-size Scaling and Universality above the Upper Critical Dimensionality 
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According to renormalization theory, Ising systems above their upper critical dimensionality 
d u = 4 have classical critical behavior and the ratio of magnetization moments Q = (m 2 ) 2 /(m 4 ) 
has the universal value 0.456947 ■ ■ •. However, Monte Carlo simulations of d = 5 Ising models have 
been reported which yield strikingly different results, suggesting that the renormalization scenario 
is incorrect. We investigate this issue by simulation of a more general model in which d u < 4, and 
a careful analysis of the corrections to scaling. Our results are in a perfect agreement with the 
renormalization theory and provide an explanation of the discrepancy mentioned. 
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One of the most important contributions to the modern theory of critical phenomena is Wilson's renormalization 
theory (see Ref. [Q for an early review). This theory explains the existence of a so-called upper critical dimensionality 
d u . It predicts that systems with a dimensionality d > d u exhibit classical exponents and violate hyperscaling, whereas 
systems with a lower dimensionality behave nonclassically. For Ising-like systems with short-range interactions, 
\Q • d u = 4. In recent years, a controversy has arisen about the value of the "renormalized coupling constant" or "Binder 
cumulant" for d > d u . On the one hand, a renormalization calculation for hypercubic systems with periodic 
boundary conditions || predicts that the Binder cumulant assumes a universal value for d > d u . On the other hand 
Monte Carlo simulations of the five-dimensional Ising model yielded significantly different results. Since the 

. renormalization theory forms the basis of our present-day understanding of phase transitions and critical phenomena, 
it is of fundamental interest to examine any discrepancies and inconsistencies with this theory. Furthermore, there 
exist several models with a lower value of d u where the above-mentioned issue may be of experimental interest 
as well. 

In this Letter, we answer the question concerning the value of the Binder cumulant. One of the key issues is the shift 
of the "critical temperature" in finite systems. We rederive this shift, which was already calculated in Ref. from 
basic renormalization equations and show that the result agrees with the shifts observed in Refs. Furthermore, 
we determine the Binder cumulant in the context of a more general Ising-like model with algebraically decaying 
interactions. This model is subject to the same renormalization equations as the aforementioned d = 5 Ising model, 
and effectively reduces to the nearest-neighbor model when the interactions decay fast enough. For slow decay, the 
upper critical dimensionality decreases below 4 and we have thus been able to investigate the question concerning the 
universality of the Binder cumulant in the classical region by means of Monte Carlo simulations of low-dimensional 
. models. This enabled us to examine a much larger range of system sizes than in the five-dimensional case. High 
statistical accuracies were obtained by using a novel Monte Carlo algorithm for systems with long-range interactions 
and we could resolve various corrections to scaling that are present. The results turn out to be in complete agreement 
with the renormalization predictions. 

We formulate our analysis in terms of the dimensionless amplitude ratio Q = (m 2 ) 2 /(m 4 ), where m is the magne- 
tization. This ratio is related to the fourth-order cumulant introduced by Binder In Ref. it is predicted that 
in hypercubic short-range Ising-like systems with periodic boundary conditions and d > 4, this quantity takes at the 
critical temperature T c the universal value 87r 2 /r 4 (i) = 0.456947 • • which is simply the value of Q in the mean-field 
model Bj. In contrast, the Monte Carlo simulations in Refs. yield the values Q w 0.50 and Q = 0.489 (6), 
respectively. In Ref. || , this discrepancy is explained by a size-dependent shift of the "effective critical temperature" 
T C (L) (defined by, e.g., the maximum in the specific heat) 

T c {L)=T c -AL- d ' 2 , (1) 

which was obtained in Refs. ^J,|| from scaling arguments. L denotes the linear system size. 

In order to examine this issue we will first outline the theoretical framework for scaling above d u . As was shown 
by Brezin J^], conventional finite-size scaling breaks down for d > d u . This is an example of Fisher's mechanism 
of dangerous irrelevant variables (see, e.g., Refs. [ [H)||it|1 ) . To examine the consequences of this mechanism for the 
finite-size scaling behavior, we briefly review the renormalization transformation for Ising-like models. Near criticality, 
one can represent the Hamiltonian for these models by one of the Landau-Ginzburg- Wilson type, 

H(^)/k B T = d d x | i(V0) 2 -h>)>+ + W 4 } • (2) 
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h is the magnetic field, ro is a temperature-like parameter and the term proportional to u keeps finite when ro < 0. 
Under a spatial rescaling with a factor b — e l the renormalization equations are, to first order in ro and it, given in 
differential form by (see, e.g., Ref. [Q) 

—jj- = Vtro + au (3a) 
^jj = , (3b) 

in which y t and y\ are the renormalization exponents of the temperature field and the irrelevant field u, respectively, 
and a is a constant depending on the dimensionality d. Upon integration, these equations yield, to first order in u, 

r' Q (b) = b yt [(r ~ cm) + aub m ~ Vt ] (4a) 
u\b) = b v 'u, (4b) 

where a is a constant. This shows that the reduced temperature t = (T — T c )/T c is proportional to (ro — au). 
Correspondingly, the free energy density / scales as 



/ It, h, u, ^) = b- d f (b yt [t + aub yi ~ Vt ], b Vb h, b Vi u, j \ + g , (5) 

where we have included a finite-size field L^ 1 and g denotes the analytic part of the transformation. The first term 
on the RHS can be abbreviated as b~ d f(t', h', u', b/L). For d > 4, the critical behavior is determined by the Gaussian 
fixed point (t*,u*) = (0,0). However, for T < T c , the free energy is singular at u = 0. Hence u is a dangerous 
irrelevant variable. The finite-size scaling properties of thermodynamic quantities can be obtained by renormalizing 
the system to size 1, i.e., setting b = L. The number of degrees of freedom then reduces to 1 and the free energy to 
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The substitution $ — (f)fv! 1 ' 4: leads to 

f(t',ti,u',l) = f(t,h), (7) 

with t = t' '/V 1 / 2 and h — h' /u' 1 / 4 . Upon renormalization, the analytic part g of the transformation also contributes 
to the singular dependence of the free energy on t, see, e.g., Ref. |12| . We absorb this contribution in the function /. 
Setting b = L and combining Eqs. (||) and yields 



/(*,M,£) =L- d f(^L"i^ 2 [t + auLy^],Ly^^) . 

For d > 4, yt = 2, yh = 1 + d/2, and y\ = 4 — d. The first argument on the RHS is the scaled temperature: 

j = L d/2 1 U + auL 2-d\ ( 9 ) 

V u 

Interpreting the term auL 2 ~ d as a shift in the "effective critical temperature" for a finite system, we recover the result 
of Ref. §. 

Let us now use the above derivation to examine the shift and rounding of critical singularities in finite systems. 
Obscrvables can be calculated from the free energy by differentiating with respect to a suitable parameter. Ignoring 
the analytic part of the free energy, we can express the thermodynamic quantities in terms of universal functions of 
the two arguments that appear in the RHS of Eq. (||). E.g., the specific heat can be written as the product of a power 
of the rescaling factor and a universal function of the scaled fields. Let the maximum of this function occur at t = c 
(c a constant). Then, the specific heat maximum occurs at a temperature which differs, in leading orders of L, from 
the critical temperature by 

At = c^L- d ' 2 ~ auL 2 - d . (10) 

The leading L-dependence of Eq. ( |To| ) agrees with Eq. (|l|). However, on the basis of Eq. (|l|) it is argued in 
Refs. 0^ that the term between brackets in Eq. (||) could be replaced by t + aL~ d / 2 , where a is a nonuniversal 



2 



constant. If this argument were correct, it would have serious consequences for the renormalization scenario: there 
must be a contribution of a new type between the square brackets in Eq. (ffa|), proportional to b~ d l 2 . There is no 
renormalization mechanism known to us which would yield such a term. Furthermore, in leading orders of L Eq. (|^) 
must be replaced by 

I = L d/2_L ( t + aL -d/A oc L d/H + a f (11) 



and in general critical-point values of finite-size scaling functions become dependent on a: they are no longer universal. 
We illustrate this for the ratio Q. Since the magnetization moments can be expressed in derivatives of the free energy 
with respect to the magnetic field, the renormalization theory predicts 

Ql(T) = Q (tV*) + qi L d -^ +■■■ . (12) 



Here Q is a universal function, t stands for the argument between brackets in Eq. (J9f), and we have introduced the 
exponents y* = yt — y\/2 and = y& — ?/i/4. The additional term q\L d ~ 2Vb = q\L~^ 2 arises from the analytic part 
of the free energy. Now suppose that Eq. ([ll]) is correct instead of Eq. (||) . Then the argument of Q is nonuniversal 
at the critical point and so is Q = lim£_,oo Ql(T c ). The value calculated in Ref. || is then just the particular value 
of Q for the mean-field model. 

Can we reconcile the renormalization scenario with the Monte Carlo results obtained up till now? The evidence for 
an effective critical temperature as in Eq. (|l|) is based upon the locations of the maxima in the susceptibility and the 
specific heat, and those of the inflection points of the absolute magnetization and the renormalized coupling constant 
cjl = —3 + 1/Ql- However, we have seen above that Eq. (||) is fully compatible with a deviation At ~ L~ d l 2 (see 
Eq. (^)). Therefore, the observed shifts do not provide evidence for the term proportional to a in Eq. ([ll]) and we 
look for a different source of the discrepancy between the renormalization and Monte Carlo results for Q. Equation (@) 



shows that there are several corrections to scaling which may well account for this. When Eq. (12) is expanded in 
tL y t , the term proportional to a yields a term qiL 2 ~ d l 2 . Furthermore, when we include a nonlinear contribution in 
u in (||), factors u in Eq. (||) are replaced by u(l + r yuL Vi ) and we find an additional term q^L 4 ^ d . Higher powers 
of these corrections may also be taken into account in the analysis, as well as the term q x L- d l 2 in @. However, 
the determination of these corrections would require accurate data for a large range of system sizes L, and the high 
dimensionality of the d — 5 Ising model presents here a major obstacle. The results presented in Refs. were based 
on 3 < L < 7 and therefore the results were by no means conclusive. Reference || used the range 5 < L < 17. Given 
these limited ranges of system sizes, it seems uncertain whether all important corrections have been resolved. Thus 
the Monte Carlo evidence against the renormalization result of Ref. || is not compelling. 

Here we follow a different approach to test the renormalization predictions. In Ref. |t]], Fisher, Ma and Nickel 
investigated the renormalization behavior of O(n) models with ferromagnetic long-range interactions decaying as 
r -(d+cr) ^ a > o). The Fourier transform of the Landau-Ginzburg- Wilson Hamiltonian is quite similar to that of 
Eq. (||); only the term proportional to k 2 is replaced by a term proportional to k a . Thus the renormalization 
equations have the same form as for short-range interactions; only the exponents and the coefficient a in Eq. (j^) 
assume different values. For < a < d/2 (d < 4), the Gaussian fixed point is stable and the critical exponents have 
fixed, classical, values (and hence hyperscaling is violated). The upper critical dimensionality is thus d u = la. In 
Fig. |l|, the regions of classical and nonclassical behavior are shown as a function of d and a. Introducing a parameter 
e = 2(7 — d, we note that the classical exponents apply for e < 0, just as in the short-range case, where e = 4 — d. In 
the limit a [ 0, each spin interacts equally with every other spin, so that we can identify this case with the mean- field 
model. Thus there is an analogy between the (short-range) Ising model with 4 < d < oo and the long-range model 
with < a < d/2. If the amplitude ratio Q has a nonuniversal value, we may therefore expect that this manifests 
itself in the long-range case as well. 

In general, the study of models with long-range interactions is notoriously difficult, due to the large number of 
interactions that have to be taken into account. However, a novel Monte Carlo algorithm § of the Wolff cluster 
type |]l3f is available that suppresses critical slowing down and, in spite of the fact that each spin interacts with every 
other spin, consumes a time per spin independent of the system size. Thus we could simulate models with algebraically 
decaying interactions in 1, 2, and 3 dimensions and obtain accuracies that were not feasible up till now (cf. Ref. |l4|| 
and references therein). For d = 1, the interaction was taken exactly Kr~( d+<1 \ whereas for d = 2 and d = 3, the 
interaction was slightly modified with irrelevant contributions decaying as higher powers of r _1 To account for 
the periodic boundary conditions, the actual spin-spin couplings are equal to the sum over all periodic images. We 
have studied linear system sizes 10 < L < 150000 for d = 1, 4 < L < 240 for d = 2, and 4 < L < 64 for d = 3, 
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generating between one and four million Wolff clusters per simulation. The ranges of system sizes are larger than in 
Refs. and more intermediate values of L are available. These facts, as well as the high statistical accuracy of 

the Monte Carlo results, allowed us to resolve the leading finite-size corrections in the Ql- 

The finite-size scaling analysis was based on the Taylor expansion of the renormalization prediction for Ql near 
criticality: 

Ql(T) = Q+pitL"* + P2 i 2 L^ + p 3 iW +■■■+ Ql L d - 2 y* +■■■ + q 3 L* + ■■■ . (13) 

The coefficients pi and qi are nonuniversal and the renormalization exponents are yt = a, y^ = (d + a)/2, and 
t/i = 2a — d. The corresponding values y£ = d/2 and y^ — 3d/ A coincide with those in the short-range case. 
In addition to the corrections to scaling in Eq. ( |l3| ) we have also included higher powers of q?,L m , which become 
important especially when a is close to d/2. In fact, omitting these corrections yielded estimates for Q close to those 
obtained in Refs. |5|||, although the residuals strongly indicated the presence of additional corrections. This confirms 
the assumption that the discrepancy between the d — 5 Monte Carlo results and the renormalization calculation is 
caused by corrections to scaling. Furthermore, the coefficient a in Eq. ([)[) is very small in all cases, in accordance 
with the fact that this correction term could not be resolved in Ref. ||. An extensive analysis of the data will be 
presented elsewhere. We have fixed all exponents at the theoretical values, in order to minimize the uncertainty in 
Q. The results presented in Table || show that the agreement between the renormalization prediction for Q and the 
Monte Carlo data is excellent. 

It could, for the purpose of comparison, be of some interest to make a correspondence between systems with short- 
range interactions in d > 4 dimensions and ci'-dimensional systems with long-range interactions decaying as r~( d +IJ \ 
Such a correspondence is possible by expressing the various finite-size scaling relations in terms of the number of 
particles N instead of the linear system size L. Then the dependence of the thermal and magnetic exponents on the 
dimensionality is adsorbed in the parameter N = L d (or L d ) and the renormalization predictions for both models 
differ only in the (modified) irrelevant exponents, (4 — d)/d and (2a — d')/d', respectively. For both models, these 
exponents vary between and —1 in the classical range, and the matching condition appears as J? = §. Hence, we 
may compare the d = 5 (short-range) Ising model with the a = |d' long-range model, i.e., a = 0.4, 0.8, and 1.2 for 
d' = 1, 2, and 3, respectively. In this sense the present work approaches the nonclassical regime even closer than 
Refs. |-|. 

Finally, we remark that models with long-range interactions provide an effective way to explore scaling properties 
above the upper critical dimensionality. E.g., the approach adopted in this Letter may be generalized to planar, 
Heisenberg and g-state Potts models, including percolation problems. For a < 2, d u is reduced by a factor a/2 in the 
case of long-range interactions. 
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FIG. 1. Dimensionality versus decay parameter a for various models. Short-range models are_described by a — 2. The open 
circles indicate the models investigated in this article and the black circle marks that of Refs. 
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TABLE I. The ratio Q and critical coupling K c for systems with long-range interactions in 1, 2, and 3 dimensions, for several 



values of the parameter a in the range < a < d/2. The numbers between parentheses represent the errors in the last decimal 
places. 
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Classical behavior 
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Nonclassical behavior 
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